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We introduce a tight-binding chain with a single impurity to act as a quantum data bus for perfect 
quantum state transfer. Our proposal is based on the weak coupling limit of the two outermost 
quantum dots to the data bus. First show that the data bus has an energy gap between the ground 
and first-excited states in the single-particle case induced by the impurity in the single particle case. 

By connecting two quantum dots to two sites of the data bus, the system can accomplish a high- 
fidelity and long-distance quantum state transfer. Numerical simulations were performed for a finite 
system; the results show that the numerical and analytical results of the effective coupling strength 
agree well with each other. Moreover, we study the robustness of this quantum communication 
protocol in the presence of disorder in the couplings between the nearest-neighbor quantum dots. 

We find that the gap of the system plays an important role in robust quantum state transfer. 
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I. INTRODUCTION 

The transfer of quantum states from one quantum unit of a quantum computer to another is of fundamental 
importance in quantum information science. Recently, in view of the great potential of a physical realization of the 
quantum computer, attention is being paid to the problem of the transfer of quantum information in a solid-state 
system. 

Recently, spin systems have been proposed as a quantum data bus for transferring information. In a pioneering 
study [l|, Bose showed that the simplest coupled spin chain with uniform nearest-neighbor (NN) couplings is able 
to act as a quantum channel, i.e., the spin system allows the transmission of an arbitrary quantum state with 
high fidelity from one end to the other. The advantages of this protocol are that no external control is required 
throughout the entire transfer process, the quantum state transfer (QST) is equivalent to the natural dynamical 
evolution of the time-independent Hamiltonian, and the system can be isolated from the environment to minimize 
decoherence. However, the drawback of this proposal is that the transfer quality decreases with the size of chain. 
One way to overcome this problem is to precisely modulate the couplings between NN spins throughout the quantum 
data bus, as suggested in Ref. 0 so as to obtain perfect QST, which is independent of the chain length. This is 
possible because the eigenvalues of the system match the parity of the corresponding eigenstates, which is a sufficient 
condition for perfect QST il. However, such an implementation requires precise control of the system, which is not 
desirable in an experiment. Another approach to achieving perfect QST is based on a gap quantum system (bl-l^. 
By weakly connecting the transmitting and receiving qubits to a gap system, the total system’s dynamics can be 
reduced to those of an effective two- or three-level system. In addition to the fact that no extra controls are required 
for communication, a key advantage of these methods is their robustness against parameter disorder, which comes 
from inevitable technological errors in the experimental implementation. Moreover, we notice that the systems with 
long-range inter-qubit interactions for perfect QST or creating entanglement are well developed as well [2a - [^ . 

In this paper, we introduce an impurity-induced gapped system (IGS), which is a tight-binding chain with on-site 
energy applied on a single quantum dot (QD), to act as a quantum channel. We demonstrate the existence of a 
nonvanishing energy gap between the ground and first-excited states in the single-particle case. We also investigate 
the QST using the IGS. It is found that at lower temperatures, the total Hamiltonian can be mapped to a three-level 
effective Hamiltonian whose energy levels are equally spaced and can be used to perform near-perfect QST. In the 
weak-coupling limit, the coupling constant of the effective Hamiltonian has an inverse relationship with the transfer 
distance. Moreover, we study the robustness of the state transfer against the static imperfections of the couplings, as 
discussed in Ref. [^, . The resulting distribution of the transfer fidelities reveals that chains with boundary states 

are more resilient to imperfections. This is reflected in more instances of high-fidelity transfer through the spin chain. 
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FIG. 1: (Color online) (a) Schematic illustrations of impurity-induced gap system (IGS) which is a A^-site chain with on-site 
energy applied on the central quantum dot. (b) A schematic setup for QST between two QDs L and R over the distance 
d — 21 + 1, throngh the gap system. 


Compared with previously proposed schemes, the advantage of our scheme is that it is simple and can be readily 
applied to experiments. 

This paper is organized as follows: In Section II, the model IGS is set up and its spectrum is introduced. In Section 
III, our QST protocol is set up and the effective Hamiltonian, Hgs, is deduced using perturbation theory. The scheme 
for using the IGS to transfer a quantum state is discussed in Section IV. Finally, conclusions of these investigations 
are presented in Section V. 


II. MODEL OF QUANTUM COMMUNICATION 
A. Quantum data bus 

We begin by introducing a one-dimensional tight-binding chain of N QDs with one diagonal impurity at Ag-th site, 
which acts as a quantum data bus. The model is shown in Fig. 1(a), which is described by the Hamiltonian 

AT-l 

Hm=-J'^ X! + + h-C.) - /To XI (1) 

3 = 1 (T=t.i cr=t,4- 

where — J (< 0) is the hopping amplitude between NN sites j and j -|-1, ^ and are the creation and annihilation 

operators of electrons on the j-th site with spin cr, and —fio (< 0) is the on-site energy of the defect. With a view 
toward the quantum information, we can encode the qubit on the spin state. Note that Eq. (1) does not contain any 
spin-spin interaction term; thus, the spin degree does not change during the evolution of the system. Hereafter, we 
shall omit the a index, denoting the electron operator with generic spin state as aj = cos 9cj^ + sin0c| This 
system can be regarded as a spinless fermion system, and the feasibly obtained results can be applied to the original 
system. In this sense, we can concentrate on the spinless fermion model in the following discussion. 

N-l 

Hm = X “ Rqo^No^No, (2) 

i=i 

For the sake of clarity and simplicity, we only consider the case where the defect is placed in the middle of the medium, 
i.e., Nq = {N + l)/2. Note that the Hamiltonian, Hm, commutes with the total number operator, h = 
and so the Hilbert space can be decomposed into subspaces corresponding to different particle numbers, n. For the 
case of transferring a single particle, we restrict the discussion to the single-particle subspace, which is spanned by 
the Fock states |j) = aj |0), with j = 1, 2,..., N. 

In this study, we focus on the bound state (or the ground state of Hm) of this Hamiltonian for nonzero fj,o, which 
can be obtained via the Bethe ansatz method. We will also show that for Hamiltonian Hm, there exists a finite energy 
gap, A. = El — Eg ^ Ro/‘^J, between the ground state and the first excited state. 

To deduce the above conclusion, we write the state in the single-particle Hilbert space as |A„) = 

Substituting the discrete superposition state into the eigenequation Hm |A„) = e„ |A„), we get 

N 

— J X fr = {R0^3,No + £n) /", 

i=l 


( 3 ) 






with open boundary condition = 0 . 

For fiQ — 0, the solution of Eq. is 
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(n + l)7rj 
2No 


( 4 ) 


with n = 0, 1, 2,..., iV — 1, the eigenvalues are = —2Jcos [(n + l)7r/2A^o]- 

We now study the effect of the impurity on the energy spectrum of Hamiltonian Hm for nonzero ^o- Before making 
calculations, we make the following observations: first, when the Hamiltonian, Hm, is processing mirror symmetry 
with respect to the inversion center, Nq, its eigenvectors, |A„), have definite parities. Moreover, if the eigenvalues, A„, 
are in increasing order, then the eigenvectors, |A„), change parity alternatively, i.e., the mirror-inverted eigenstates, 
|A„), satisfy the relation /" = upon assuming that even (odd) n label even (odd) eigenstates |A„). 

Second, the probability density of all the eigenstates with odd parity in the central site, Nq, is zero, i.e., = 0, 

which means that the eigenstates with odd parity are unaffected by the presence of the impurity. Third, by the 
Hellmann-Feynman theorem, the eigenvalues of even-parity eigenstates decrease due to the presence of the impurity. 
Furthermore, the impurity contributes exactly one bound state, which we focus on in this study. 

To see more precisely what happens for /tq 7 ^ 0, we solve Eq. ([3]) via the Bethe Ansatz method. In this study, the 
bound state is the ground state of Hm ■ Through a straightforward calculation, one can obtain the following analytical 
result for the ground state 




sinhfcoj, j < Nq 

sinhfco(A^ + 1 - j), j > Nq 


( 5 ) 


with the eigenvalue Aq = + 1 , where kQ = 

g2koNo ^g2feo y'4 ^g 2 feo _ x) is the renormalization factor. 


In 


? + +1 and ^ = /ro/2J; Q 


The remaining eigenstates with even parity are extended and similar to Eq. O; the appropriate Ansatz is 


which yields the eigenvalue e„ 


f sin knj, j < Nq 

\ sin kn{N +1- j), j > Nq ’ 


—2JcosA:„ and the wave vector, kn, obeys 


^ sin knNo = cos knNo sin fc„ 


( 6 ) 

( 7 ) 


Setting tan</ 5 „ = sinfc„, Eq. ([7|) becomes cos {knNo + (pn) = 0, whose allowed values are 


kji 


(2m — 1 ) TT — 2ipn 
2An 


,m = 1,2..., A^o - 1- 


( 8 ) 


From the above equations, we know that (i) the phase shift = 0 for ^ = 0 and = tt/2 for ^ = 00 , and that (ii) 
the phase shifts do not alter the order of the sequence {kn}. 

Until now, we have only discussed the solutions of eigenequation Hm |A„) = £„ |A„) without any external pertur¬ 
bation. In the thermodynamic limit where Nq — >■ 00 , the excited energies become a continuous energy band; it is not 
hard to hnd that the energy gap between the ground state and the first excited state (see the Fig. [3^) is 


A = 2 Jv '^2 + 1 _ 2J. 

For very small values of onsite energy, i.e., fiQ <C J, we get A Ri 


(9) 


B. The subspace Hamiltonian 

Now, let us introduce the protocol of quantum communication by using IGS, in which two individual QDs (sender 
and receiver) are symmetrically coupled to an IGS on opposite sides of the data bus (as illustrated in Fig. 1(b)). 
Moreover, QDs L and R are supplied with on-site energy, —p. The total Hamiltonian consisting of (N + 2) QDs reads 

H = 

Ho = 

H' = 


Hq + iJ', 

Hm - Ai [a'^L^L + 


-Jo fa^dATo-i + d^djvo-i-i + h.c.j 


) 


( 10 ) 
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FIG. 2: Schematic illustration of the energy levels of the system, (a) When the connections between two QDs and the medium 
switch off Jo = 0 the ground states are three-fold degenerate, (b) When Jo = 0 switches on, the degeneracy will be removed 
and split into three sub-levels with the level spacing S = Jeft. This is approximately equivalent to that of tight-binding chain 
with three QDs. 


where ol and ur are the annihilation operators of electron on L and R, {Nq ± 1) denotes the connecting sites of the 
chain, and the coupling constant, Jq, measures the strength of the interaction. 

In the absence of coupling between the two qubits and the medium (Jo = 0) and setting —/i = Ao, the total 
Hamiltonian (ITOl) can be diagonalized in the basis {| J), \R) , |Ao), |Ai) ,..., IAat-i)}, and its ground states are threefold 
degenerate, i.e., |L), |Ao), and \R) have the energy Eg^'^ = — /r. These three states can be regarded as the components 
of an effective pseudo-spin-1 system that span an invariant subspace. The original energy degeneracy will break down 
by switching on the weak coupling, Jq {Jq <C A), and the ground state will split into three sub-levels with level 
spacing AE = | Jeff|, as illustrated in Fig. HKb). Here, Jeff is the effective hopping integral that can be calculated 
as follows. 

When switching on Jq, the eigenequation becomes {Hq + H') |'0) = E For weak coupling between QDs and 
the bus, Jo <C A, El' can be treated as a perturbation Hamiltonian. Let us assume that, in some definite way, 
we can divide the basis into two classes, [G] = {|L), |i?), |Ao)} where the relative complement of [G] is denoted by 
[O] = {|Ai),..., IAat-i)}. Defining 

g = |J)(J| + |Ao)(Ao| + |i?)(i?|, 

N-l 

n— 1 

denote two orthogonal projection operators of two different subspaces. It is easy to check that QO = 0 and satisfying 
g + O = I. The eigenequation can be rewritten as 

{g + d'^H(g + 6 ) (g + o) \^p) = E(g + d'j \^). (13) 

The above equation can be decomposed into two basic formulae in subspaces (G) and (O) 

Hggg W + HgoO = Eg\i;), 

Hogg IV’) + HooO \ip) = EO \ii )), 

where Hap = aHjl, ^q;,/ 3 = ^,(0^ Using Eq. (HU, one can express G |'0) in terms of ^ |'0): 

d\il;)=(E-Hoo)~^ Hogg\ilj). (16) 

so that, substituting the above equation into Eq. da, one finds that, to second order, the equation only evolves 

QW- 


(14) 

(15) 


( 11 ) 

( 12 ) 


H,sg H) = Eg 1 ^), 


( 17 ) 
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where 

HeS = Hg-p + Hpo (yE — Hoo'j Hog (18) 

denotes the effectvie Hamiltonian in subspace (A) with 

Hgg = -JoCo{\L) + \R)){Xo\+h.c.-^ig, (19) 

N-l 

Hoo = ^A„|A„)(A„|, (20) 

n—1 

N-l 

Hgo = ^-JoCn(0[|i) + (-in^)](A„|. (21) 

n—1 

and Cn (0 = {Hq — ^ |A„). Through a straightforward calculation, one can obtain 


Hgo (^fJ- + Hoo^ Hog = 


N-l 

E 

n—1 


JS\U 

E-Xr 


[|L) + (-l)”|i?)][(L| + (-l)"(i?|]. 


Note that the eigenvalues, E, determined from Eq. (EZD, are perturbed eigenvalues around respective unperturbed 
value With this connection, one seldom requires the second-order correction, which is small (Jq <C I-E — A„|, 
which is the condition for the perturbation procedure to be a good approximation in this problem); it is therefore 
sufficient to quote the first-order results 


Eeff - Jeff (|E) + m (Aol - I (|L) (L| + |Ao) (AoI + \R) {R\) + h.c., (22) 


with effective coupling strength Jeff = JqCo (0- 

In this section, we have shown that the total Hamiltonian (HOI) can be simplified to the effective Hamiltonian (I22p . 
due to a large gap (compared with coupling strength Jq) existing in the medium. This approximation holds when the 
energy splitting, v^Jeff, caused by the iJeff is smaller than the typical gap for the unperturbed Hamiltonian, Eq; i-e., 
Jeff <C A. To check the range of validity of the above effective Hamiltonian, we compare the analytical result of Jeff 
with the results {Ei — Eq) /^/2 obtained by direct numerical diagonalization of the Hamiltonian (fTOll . The results of 
this comparison are plotted in Fig. 3 for a system of iV = 499, with Jg = 0.001 J, and = 0.1 J, 0.05J, and 0.01 J. 
In this figure, one can see that taking ^o/Jo bigger than 50, the effective coupling strength, Jeff, of Heff agrees very 
well with that obtained numerically. So far, the validity of the effective Hamiltonian (1^^ is firmly established. Thus 
one should be able to obtain high-fidelity QST with the effective Hamiltonian whenever the perturbation solution is 
valid. Furthermore, we will show that the existence of an energy gap can also be used to protect the performance of 
QST in the presence of static disorder in the couplings of the quantum data bus. 

However, it is worth pointing out that large /ig can improve the validity of iJeff but decrease the transfer efficiency 
characterized by Jeff, since l/JeS determines the transfer time of the QST between the two qubits, L and R. As 
observed in Fig. 3, the decay rate of Jeff directly depends on the value of ^g. The smaller the /rg is, the slower the 
decay rate will be. Typically, Jeff decreases almost linearly with the increase of the transfer distance for /rg = 0.01 J. 
From the two competing aspects described above, we can summarize the proper choice of the system parameters, ^g 
and Jg, for high-fidelity QST. 

To briefly summarize, we have theoretically and numerically studied Jeff as a function of d in a specific range 
of parameters. However, the obtained conclusion is based on the fact that the id^gff given by Eq. (1221) is a valid 
approximation in the studied range. In the following discussion, the validity of iJeff is investigated by comparing the 
eigenstates of HeS with the lowest three states of the total system (uni). 

Define the quasi-angular momentum states |j, m) as 

|1,0) = ^{\L)-\R)), (23) 

|1,±1) = i(|L)±^|Ag) + |i?)), (24) 

which are the eigenstates of effective Hamiltonian ([2^ . On the other hand, the eigenstates of H can be generally 
written as 


N-l 

|E) ~\- ^ ^ Cji [Att,) -|- Cjl |-/?) , 
n—D 


(25) 
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FIG. 3: (Color online) (a) Behavior of JbS as a function of the transfer distance d for a chain of A'’ = 499 sites, Jo = 0.001 
and different values of Curves from top to bottom are, respectively, for fio = O.IJ, 0.05J, and O.OIJ. Continuous curves 
display the approximate expression for the Jeff and symbols denote the exact numerical data which is given by {Ei — Eg) ly/2. 
(b) The same as in (a) but for /ro = 0.05J and different values of Jo. 


where we have the condition |cl|^ + |cn|^ + |cfi|^ = 1 for the normalization of |'0jm)- Moreover, we assign the state 

to denote the ground state for j = 1, m = 1, the first excited state for j = 1, m = 0, and the second excited 
state for J = 1, m = — 1. To evaluate the fidelity of the Jfeff induced by the perturbation, we introduce the overlap 

Pjm = IUjITT- • (26) 

For the case where Jq = 0, the ground states \ijjjm) of H are threefold degenerate and can be written in 

symmetrical form by linear combinations of |L), |Ao), and \R). Under this condition, one can obtain Pjm = 1 for j = 1 
and TO = 0, ±1. In particular, we have |cl|^ = |c_r|^ = 1/4, \cof = 1/2 for to = ±1 and \cl\^ = \cr\^ = 1/2, |co|^ = 0 
for TO = 0. For the practical Hamiltonian JJ, i.e., Jq ^ 0, the values of |ci|^ (i = L,0,R) and Pjm are numerically 
calculated for the three lowest eigenstates in the N = 499 system with /ip = O.IJ, 0.05J and Jq = 2 x 10“^J for finite 
transfer distances d = 5, , 15, 25, 35, 45, 55, and 65, which are listed in Tables 1(a) and (b). 

We remark that the condition for mapping H to the effective Hamiltonian (1221) is that Jq must be small enough 
compared to the energy gap A of the medium rather than the on-site energy fiQ. As mentioned before, the energy 
gap is A Ri Mo/dJ for small /tq (compared with J). It is straightforward to obtain A r; 2.5 x 10“^J for = O.IJ 
and A r 6.25 x 10“^ J for fiQ = 0.05 J. From the numerical results shown in Table I, we observe that the realistic 
interaction leads to the results for |ci|^ {i = L, 0, R), which are very close to those described by iJeff, even if A is of the 
same order of Jq. It is clear that such a three-level subsystem allows state \L) to transfer with high fidelity, and the 
coherent population exhibits oscillations between the QDs on the two ends. The oscillation period of the population 
is given by r = 7r/-\/2Jeff, and we can say that the quantum state is transferred from QD L to QD R at the time 
T = (2n — 1) X r. 


III. QUANTUM STATE TRANSFER 
A. Weak Coupling Regime 

Note that the spectrum structure and the corresponding parity of the effective Hamiltonian, iJeff, obey the spectrum- 
parity matching condition (SPMC) Si exactly, which is the general criterion for perfect QST. In this section, we 
consider the QST scheme based on our system. Assume Alice is at the sender site. A, and Bob is at the receiver 
site, B. Let Alice hold an electron with a spin state that she wants to communicate to Bob of \<f) = cos (0/2) |t) + 
sin (0/2) II), where |t) (||)) denotes the spin-up (down) state. Thus, the initial state of the total system is 

|4> (0)) = |L) = ^cos0c^ .J. -I- e®'^sin0c^ |0)^, which is a superposition of the eigenstates of Hamiltonian 

l'f'(0)) = i(|l,l) + V2|l,0) + |l,-l)). 


(27) 
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TABLE I: The overlap Pjm and its three components, which provide a criteria for the validity of are calcnlated nnmerically 
for the gronnd state, first excited state and second excited state of total system for finite transfer distance d = 5, 15, 25, 35, 
45, 55, and 65. The results for /ro = O.IJ, and 0.05J (Jo = 2 x 10“®J) are listed in (a), and (b) respectively. It shows that the 
result based on the realistic interaction is very close to that by H^s even if /io is not large enough. 


States 

j 

m 


5 

15 

25 

35 

45 

55 

65 

(a) 



|ci.r 

0.2552 

0.2531 

0.2539 

0.2569 

0.2618 

0.2687 

0.2778 

I'Aii) 

1 

1 


0.4884 

0.4932 

0.4920 

0.4860 

0.4757 

0.4613 

0.4426 





0.2552 

0.2531 

0.2539 

0.2569 

0.2618 

0.2687 

0.2778 




Pll 

0.9986 

0.9994 

0.9997 

0.9995 

0.9988 

0.9973 

0.9950 





0.4999 

0.4994 

0.4987 

0.4980 

0.4975 

0.4971 

0.4968 

I'Aio) 

1 

0 

Icol^ 

3.457x10"^® 

2.334xl0"^'‘ 

6.501x10“^® 

2.248x10“^® 

1.044x10“^® 

2.920x10"^® 

1.992x10"^^ 





0.4999 

0.4994 

0.4987 

0.4980 

0.4975 

0.4971 

0.4968 




Pio 

0.9999 

0.9988 

0.9973 

0.9960 

0.9949 

0.9941 

0.9935 





0.2432 

0.2462 

0.2459 

0.2429 

0.2375 

0.2301 

0.2205 


1 

-1 

Icol^ 

0.5116 

0.5068 

0.5080 

0.5140 

0.5243 

0.5387 

0.5574 





0.2432 

0.2462 

0.2459 

0.2429 

0.2375 

0.2301 

0.2205 




pl-1 

0.9979 

0.9992 

0.9997 

0.9995 

0.9987 

0.9973 

0.9950 

(b) 




0.2641 

0.2599 

0.2580 

0.2580 

0.2595 

0.2625 

0.2666 

li’ii) 

1 

1 

Icol^ 

0.4633 

0.4741 

0.4801 

0.4817 

0.4795 

0.4738 

0.4650 




IckI^ 

0.2641 

0.2599 

0.2580 

0.2580 

0.2595 

0.2625 

0.2666 




Pll 

0.9904 

0.9933 

0.9957 

0.9973 

0.9981 

0.9980 

0.9970 





0.4999 

0.4985 

0.4963 

0.4937 

0.4912 

0.4887 

0.4865 

I'Aio) 

1 

0 

Icol^ 

6.467xl0"^'‘ 

2.461x10"^® 

1.581x10“^® 

2.075x10“^® 

2.362x10“^® 

7.620xl0"^'‘ 

2.793xl0"^'‘ 




IckI^ 

0.4999 

0.4985 

0.4963 

0.4937 

0.4912 

0.4887 

0.4865 




Pio 

0.9997 

0.9970 

0.9925 

0.9875 

0.9823 

0.9774 

0.9730 




Icrl^ 

0.2194 

0.2286 

0.2348 

0.2379 

0.2383 

0.2360 

0.2315 

I'Ai-i) 

1 

-1 

Icol^ 

0.5318 

0.5231 

0.5185 

0.5177 

0.5203 

0.5261 

0.5350 




IckI^ 

0.2194 

0.2286 

0.2348 

0.2379 

0.2383 

0.2360 

0.2315 




pl-1 

0.9684 

0.9791 

0.9874 

0.9931 

0.9963 

0.9974 

0.9967 


At time t, the initial state Iflf (0)) evolves into 

\-qj{t)) = IT-(0)) 

= i(e*^*|l,l) + V2|l,0)+e-*^‘|l,-l)) (28) 

where 6 = and we have neglected the overall phase, The density matrix corresponding to I'll (t)) is 

p =\^ (t)) ('ll (t)|, and the probability of state |'E' (0)) transferring to the QD R at time t is defined as 

F{t) = Tt{ppr) = sin'‘ . (29) 

At the moment when t = t = tt/S, F{t) = 1 indicates that our scheme can perform QST perfectly. That is to say, 
the system evolves into a new factorized state 

l^(r)) = i(e^-|l,l) + y2|l,0) + e—|1,-1)) 

= e^’^li?). (30) 

As an example of verifying the validity of the effective Hamiltonian HeS, the fidelity for N = 499 and transfer 
distance d = 5, with Jq = 2 x 10“^J, and po = O.IJ, 0.05J are plotted in Figs. 4(a) and (b). They show that small 
Jo leads to a result for transfer fidelity, which is very close to that described by the effective Hamiltonian, Heff- 






FIG. 4: (Color online) Comparison of the exact transition probability F{t) (open circle) with the analytic result of Eq. (27) 
(red solid line) for the system with N = 499, 1 = 2, Jo = 2 x 10~®J, /io = O.IJ (upper plot), and 0.05J (lower plot). Time is 
expressed in units of 10‘*/J. It shows that, small Jo leads to the result about transfer fidelity is very close to that described by 
the effective Hamiltonian iJeff. 


B. Robustness to Disorder 


We now turn to the performance of spin chains in the presence of static imperfections in the couplings, which are 
unavoidable in experimental implementations. We will show that the energy gap can protect the performance of the 
QST in the presence of static disorder in the system parameters. 

We now assume that the tunnel coupling of the medium Hamiltonian has a random but constant offset, 6ej, i.e.. 


H' = — J(1 + 5ej) (a]dj+i + h.c.) 

f a[aL + j - Jq \^a^j^aNo-i + a^naNo+i + h.c. j , (31) 


where <5 is the maximum coupling offset bias relative to J; ej is drawn from the standard uniform distribution in the 
interval [— 1 , 1 ] and all ej are completely uncorrelated with all sites along the chain. 

We numerically calculate the Schrodinger equation for the dynamical evolution and compute the overlap, F{t) = 
2 


(i?| e ‘ \L) , to assess the performance of the chain. In Fig. [5] we plot the behavior of F{t) as a function of time, 

t, in the system with N = 499 QDs, I = 2, Jq = 2 x 10“^ J for three cases: (a) fio = O.IJ, 5 = 5 x 10“^, (b) fio = O.IJ, 
(5 = 1 X 10“^, and (c) fiQ = 0.5J, (5 = 1 x 10“^. From this comparison, one can see that (i) this scheme is robust 
against the static disorders that would be unavoidable in experimental implementations and (ii) that the large energy 
gap (or large fio) is more robust than small one against disorder. 


IV. SUMMARY 

According to quantum mechanics, it is not difficult to establish a long-distance QST using a gap system. However, 
the magnitude of the gap in this kind of scheme is crucial: first, the gap should be independent of the size of the 
system; second, the energy gap should be manipulated as required for perfect QST. The reason is that if the gap is 
too large, the QST period increases exponentially with the distance between two distant parities; when the gap is too 
small, the fidelity of the QST is reduced. 

In this study, the quantum transmission of an electron through an IGS (serving as the data bus) is studied by 
theoretical analysis and numerical simulation. First, we show that the IGS has a nonvanishing energy gap above the 
ground state, which depends only on the on-site energy, —/io, of the impurity. The approach to realize perfect QST is 
based on weakly connecting two external QDs with the bus. Different transfer distances can be achieved by suitable 
choices of connecting sites to the data bus. By treating the weak coupling as a perturbation, we find that a gap 
system can induce an effective three-level Hamiltonian [Eq. (I22p ]. This theoretical result is confirmed by performing 
numerical simulations; moreover, the effective coupling, Jeff, also decays slowly with increasing transfer distance if the 
system parameters are chosen reasonably. 
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FIG. 5: (Color online) The transition probability F(t) of QST as a function of time in a = 499 system with Jo = 2 x 10~® J, 
I = 2 and random imperfections of coupling strength S function of time. The left figure corresponds to the case when /io = O.IJ 
and (5 = 5 X 10“®, the middle figure to /io = O.IJ and <5 = 1 x 10“^, and the right figure to /to = 0.5J and 5 = 1 x 10“^. The 
results shows that weak fluctuations in the coupling strengths do not deteriorate the quality of QST due to the exitance of 
energy gap. 


Furthermore, the fault tolerance for more realistic system parameters is also demonstrated. It has been shown that 
perfect state transfer can also be achieved in the presence of disorder. For larger values of the energy gap (or /tq), the 
effect of disorder on the quality of QST will be strongly suppressed. 
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